Molecular motors stiffen non-afRne semiflexible polymer networks 
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Reconstituted filamentous actin networks with myosin motor proteins form active gels, in which 
motor proteins generate forces that drive the network far from equilibrium. This motor activity 
can also strongly affect the network elasticity; experiments have shown a dramatic stiffening in in 
vitro networks with molecular motors. Here we study the effects of motor generated forces on the 
mechanics of simulated 2D networks of athermal stiff filaments. We show how heterogeneous internal 
motor stresses can lead to stiffening in networks that are governed by filament bending modes. The 
motors are modeled as force dipoles that cause muscle like contractions. These contractions "pull 
out" the floppy bending modes in the system, which induces a cross-over to a stiffer stretching 
dominated regime. Through this mechanism, motors can lead to a nonlinear network response, even 
when the constituent filaments are themselves purely linear. These results have implications for 
the mechanics of living cells and suggest new design principles for active biomemetic materials with 
tunable mechanical properties. 



The mechanics of living cells is largely governed by 
the cytoskeleton, a complex assembly of various filamen- 
tous proteins. Cross-linked networks of actin filaments 
form one of the major structural components of the cy- 
toskleton. However, this cytoskeleton is driven far from 
equilibrium by the action of molecular motors that can 
generate stresses within the meshwork of filaments [1-3]. 
Such motor activity plays a key role in various cellular 
functions, including morphogenesis, division and locomo- 
tion. The nonequilibrium nature of motor activity has 
been demonstrated in simplified reconstituted filamen- 
tous actin networks with myosin motors [4-8]. Even in 
the absence of motor proteins, such in vitro networks of 
cytoskelctal filaments already constitute a rich class of 
soft matter systems that exhibit unusual material prop- 
erties, including a highly nonlinear elastic response to 
external stress [9-15]. This nonlinear response can be 
exploited using molecular motors [4, 7]; the network stiff- 
ness can be varied by orders of magnitude, depending on 
motor activity. A quantitative understanding of such ac- 
tive biological matter poses a challenge for theoretical 
modeling [3, 16-21]. 

The nonlinear mechanical response of reconstituted 
biopolymer networks in many cases reflects the nonlinear 
force-extension behavior of the constituting cross-links 
or filaments [9-11, 14, 22]. For such networks, there is 
both theoretical and experimental evidence that internal 
stress generation by molecular motors can result in net- 
work stiffening in direct analogy to an externally applied 
uniform stress [4, 7, 18-20, 23]. However, the mechani- 
cal response of semiflexble polymers is highly anisotropic 
and is typically much softer to bending than to stretch- 
ing. In some cases, this renders the network deforma- 
tion highly non-affine with most of the energy stored in 
bending modes [24-28] . Such non-affinely deforming stiff 
polymer networks can also exhibit a nonlinear mechan- 
ical response, even when the network constituents have 
a linear force-extension behavior [29-32]. However, the 
effects of internal stresses generated by molecular motors 
in such networks are unknown. 




FIG. 1: Example of a portion of the diluted 2D phantom tri- 
angular network at Q = 1/4 and K = 10~ 3 . The freely hing- 
ing binary cross-links are indicated in black. Motors generate 
muscle-like contractions, which we model with force dipoles. 
The segments along which these contractile force dipoles act 
are indicated with red dumbbells. The inset shows an en- 
largement of the network. 



Here we study the effects of motor generated forces on 
the network mechanics in 2D networks of athermal, stiff 
filaments using simulations. In the absence of motors, 
these networks can exhibit strain stiffening under an ex- 
ternally applied shear. This behavior has been attributed 
to a cross-over between two mechanical regimes; at small 
strains the mechanics is governed by soft bending modes 
and a non-affine deformation field, while at larger strains 
the elastic response is governed by the stiffer stretch 
modes and an affine deformation field [29]. We show that 
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motors that generate internal stresses can also stiffen the 
network. The motors induce force dipoles leading to mus- 
cle like contractions, which "pull out" the floppy bending 
modes in the system. This induces a cross-over to a stiffcr 
stretching dominated regime. Through this mechanism, 
motors can lead to network stiffening in non-afline stiff 
polymer networks in which the constituting filaments in 
the network are themselves linear elements. These re- 
sults have implications for the mechanics of living cells 
and propose new design principles for active biomemetic 
materials with highly tunable mechanical properties. 

I. THE MODEL 

To study the basic effects of internal stress gener- 
ated by molecular motors on the macroscopic mechanical 
properties of stiff polymer networks we employ a mini- 
malistic model, which is illustrated in Fig. 1. Filamen- 
tous networks in 2D are generated by arranging filaments 
spanning the system size on a triangular lattice. Since 
physiological cross-linking proteins typically form binary 
cross-links, we randomly select two out of the three fil- 
aments at every vertex between which we form a binary 
cross-link. The remaining filament crosses this vertex as 
a phantom chain, without direct mechanical interactions 
with the other two filaments. The cross-links themselves 
hinge freely with no resistance. With this procedure we 
can generate disordered phantom networks, based on a 
triangular network, but with local 4-fold (z — 4) con- 
nectivity corresponding to binary cross-links. The use of 
a triangular lattice avoids, for example, well-known me- 
chanical pathologies of the 4-fold square lattice. To cre- 
ate quenched disorder in the network, we cut and remove 
filament segments between vertices with a probability Q. 
This also has the effect of shortening the filaments. 

The filaments in the network are described by an ex- 
tensible wormlikc chain (EWLC) model with an energy 

where k is the bending rigidity, t is the tangent vec- 
tor at a position s along the polymer backbone and 
is the local relative change in contour length, or 
longitudinal strain. We can quantify the relative im- 
portance of the stretch and bend contributions by the 
lengthscale £b = y/n/fi; this length scale forms one of 
the key control parameters for the network mechanics. 
For simple cylindrical beams with a radius r, the stretch 
modulus \i is related to k through ^ mec h = 4k/ r 2 , and 
£fe = r/2. In contrast, a thermally fluctuating semi- 
flexible polymer segment cross-linked in a network on 
a length-scale £ c also has an entropic thermal stretch 
modulus yu t h = 90K 2 /fceT^ [33], where ks is Boltz- 
mann's c onstant and T is the temperature. In this case, 
If, = £ c iJ£ c /90£ p , where £ p — n/k^T is the persistence 
length. The most relevant values of £b/£ c for biopolymer 



systems range from 10~ 2 — 10 _1 . This range extends 
from relatively stiff actin filaments to the more flexible 
intermediate filaments. Various actin binding proteins 
are capable of forming tightly coupled stiff bundles of 
actin filaments, which further reduces £b- The mechanical 
and thermal moduli add as springs in series and the total 
modulus is given by n^ 1 = /Vech+^th • ^ n ^ ne remainder 
of this paper all lengths are determined in units of the dis- 
tance between lattice vertices £q and the bending rigidity 
k is measured in units of \i£\. Here, we focus on non- 
linearities arising in networks of purely linear elements. 
Thus, we do not include intrinsic nonlinearities associ- 
ated with the force-extension curve of thermal filaments. 
This has been examined theoretically in Refs. [18-20] 

In our numerical simulations we use a discretized ver- 
sion of Eq. (1) with a node at and between every lattice 
vertex. The mid-node allows us to capture buckling down 
to the single segment length-scale. To model the effect 
of muscle like contractions induced by molecular motors, 
we introduce force dipoles in the network [4, 18, 19, 23]. 
These force dipoles are randomly placed at neighboring 
cross-links. The force dipoles only act along existing 
bonds and, therefore, do not introduce additional con- 
straints in the network. The total energy of the sys- 
tem includes a sum of the EWLC Hamiltonian over all 
filament segments and the work extracted by the force 
dipoles 

£ = E^-E/^-< ( 2 ) 

i <ij> 

where is the distance between cross-link i and j. The 
force dipoles are numerically implemented by shorten- 
ing the effective rest length of the bond along which the 
motors acts in the stretch term of the energy (Eq. 1). 

The rest length is reduced by an amount 5r^; the re- 
sulting force is given by fiSr^ /£ < H- The effects of 
internal motor generated stresses modeled in this way is 
illustrated in Fig. 1. 

To investigate the mechanical response of the network, 
an external strain 7 is applied by translating one of 
the horizontal boundaries to which the filaments are at- 
tached. The internal degrees of freedom of the network 
are relaxed by minimizing the energy using a conjugate 
gradient algorithm [34]. To reduce edge effects periodic 
boundary condition are employed at all boundaries. The 
linear shear modulus of a network of size W 2 is related 
to the energy G — ^2^2 for small strains. In the non- 
linear regime it is common to determine the differential 
modulus K = ppa^f , which reduces to G for small 7. 
Similarly, the stress can be calculated in the nonlinear 
regime through cr ext = p^^- These measurements al- 
low us quantify the mechanical response of the system. 
Here we use system sizes ranging from W 2 ~ 3000 to 
8000. 
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FIG. 2: The differential shear modulus K = Aa/A r y as a func- 
tion of the applied external stress <r cxt for various ratios of 
bending rigidities re and fixed Q = 1 /4. K and <T e xt are mea- 
sured in units of fi/io. The affine prediction is shown as a red 
dashed line, which constitutes an upper bound to the elastic 
response. Although definite powerlaw regimes appear to be 
absent, the stiffening curves for floppy systems with re < 10 -3 
initially show a stiffening behavior of approximately K ~ a 
that crosses over to a regime K ~ er 1//2 at large stress, as 
shown by the dashed lines indicating slopes of 1 and 1/2. For 
stiffer systems with k > 10~ 2 , only the second of these regimes 
is apparent. The inset shows the linear shear modulus G as 
a function of k, and the red dashed line indicates the affine 
prediction. 

II. RESULTS AND DISCUSSION 
A. Passive networks 

We probe the 2D phantom triangular networks by de- 
termining both the linear and nonlinear elastic response 
of the networks in the absence of motors. The linear me- 
chanical response of diluted networks (Q > 1) exhibits 
two distinct mechanical regimes. At low re, the shear 
modulus G scales directly with shown in the in- 

set of Fig. 2. This demonstrates that in this regime the 
macroscopic mechanics is governed by filament bending 
deformation modes. By contrast, at large re the shear 
modulus asymptotically approaches a limit in which G 
is independent of re indicative of a stretching dominated 
regime. These result are consistent with previous obser- 
vations on 2D mikado networks [24-26]. 

These mechanical regimes have important implications 
for the nonlinear elastic response. When a large external 
shear is imposed on a network that is initially in the bend- 
ing dominated regime, the differential modulus K = ^ 
increases strongly as a function of external stress <7 cxt , as 
shown in Fig. 2. Previous studies have observed similar 
stiffening in networks with strictly linear elements [29- 
32]. This remarkable behavior has been explained in 
terms of a siram-induced cross-over from a bending to 
a stretching dominated regime. At low stresses the net- 
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FIG. 3: The differential shear modulus K = dcr/d7 as a func- 
tion of the applied external stress a ex t for various values of (L) 
at fixed bending rigidity re = 10 -3 . K and <T ex t are measured 
in units of fi/£o- Although definite powerlaw regimes appear 
to be absent, the stiffening curves for (L) < 5 initially show a 
stiffening behavior of approximately K ~ a that crosses over 
to a regime K ~ a 1 ^ 2 at large shear, as shown by the dashed 
lines that indicate a slope of 1 and 1/2. For longer filaments, 
only the second, weaker stiffening response is apparent. The 
inset shows the linear shear modulus G as a function of the 
average filament length (L) , and the red dashed line indicates 
the affine prediction in the high molecular weight limit. 

work mechanics is governed by bending modes, which 
for small re constitute the softest modes in the system. 
However, when the stress is increased the deformations 
become correspondingly large and the stretching of fila- 
ments is no longer avoidable. This picture is consistent 
with our simulations. When a substantial shear is im- 
posed the stiffening curves — over a large range of bend- 
ing rigidities — converge to a single curve that is consis- 
tent with the affine prediction, shown as a red dashed 
line in Fig. 2. This calculation also demonstates that 
even an affinely deforming network of strictly linear el- 
ements stiffens under shear. This stiffening behavior is 
purely due to geometric effects; under shear the network 
becomes increasingly anisotropic and the filaments reori- 
ent to line up in the shear direction [35]. The extent of 
this purely geometric stiffening is, however, limited, as 
can be seen in the figure. Moreover, this geometrically- 
stiffened limit represents an upper bound on the stiffness 
of networks with purely linear elements. Such systems 
cannot stiffen indefinitely. 

In addition to re, the average length of filaments in 
the system (L) constitutes an important control param- 
eter for the linear response. We can probe this by vary- 
ing Q, since the average length of filaments is given by 
(L) = 1/Q [36]. Consistent with previous work [24-26], a 
cross-over from a non-affinc bending regime and an affine 
stretching regime can also be achieved by increasing (L), 
as shown in the inset of Fig. 3. In the high molecu- 
lar weight limit, (L) — > oo, the system responds purely 
affinely. We estimate that in experimental biopolymer 
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FIG. 4: The shear modulus G as a function of force exerted 
per motor fo for various motor densities pM at fixed 2 = 1/4 
and k = 10~ 3 . The shear modulus G is normalized by the 
shear modulus Go of the passive network. The inset shows 
the shear modulus Go as a function of the generated stress 
ctm- The apparent collapse of these curves supports supports 
the hypotheses that cm is the appropriate control variable. 

systems (L) varies a over a range of order 5-30, in units 
of the network mesh size. The strong dependence of the 
linear elastic response on (L) is also reflected in the non- 
linear response (Fig. 3). Networks with shorter filaments 
are increasingly governed by soft bending modes and thus 
exhibit a greater degree of stiffening under shear. 

In the absence of motors, we find that our diluted 
phantom triangular networks exhibit a linear and non- 
linear response to external shear that is consistent with 
previous work on 2D off-lattice networks of stiff fila- 
ments [24-26]. Our phantom triangular networks thus 
provide a good model system to study the effects of inter- 
nal stresses generated by molecular motors in athermal 
networks. 



B. Active networks 

To investigate the effect of motor generated stresses we 
introduce force dipoles in the network at various densi- 
ties pm- The shear modulus G increases strongly when 
the force exerted by a single motor f is increased beyond 
a threshold value, as shown in Fig. 4. Interestingly, the 
motor forces at which the system becomes nonlinear for 
low motor densities is close to the buckling force thresh- 
old fb = tt 2 k/£c Ri 2 x 10~ 3 . The buckling force threshold 
has been identified as an important force-scale for stiff- 
ening of these networks under external shear [29, 32]. In 
addition, these data imply that a minimum motor den- 
sity is required for motor generated stiffening, consistent 
with recent experiments [7]. The characteristic motor- 




FIG. 5: The linear shear modulus G as a function of motor 
generated stress <tm for various of bending rigidities re at fixed 
Q — 1/4. The stiffening curves for networks with floppy fil- 
aments for re < 10 -3 show an approximate scaling behavior 
given by K ~ a, as shown by the dashed lines that indicate a 
slope of 1. The inset shows the critical stress for the onset of 
stiffening as a function of re for both the active (red squares) 
and the passive (black circles) systems. 

generated stress can be expressed as om = Pm^o/o- Re- 
markably, all stiffening curves can be collapsed by ex- 
pressing the shear modulus as a function of ctm (upper 
inset Fig. 4). This demonstrates that the characteris- 
tic motor generated stress ctm is a useful quantity, even 
though the distribution of stress is likely to be highly 
heterogenous. 

To explore the nature of the stiffening induced by mo- 
tors we study the networks' response at various values of 
re. We observe that motor activity dramatically increases 
the network stiffness over a range of re values, as shown 
in Fig. 5. Interestingly, the degree of stiffening induced 
by motors stress is substantially larger for networks with 
lower re, while for large re we observe no stiffening at all. 
To compare the stiffening between the active and passive 
networks, we determine the critical stress for the onset 
of stiffening. When the linear mechanics of the networks 
is controlled by bending modes (G ~ re) we find that a c 
scales linearly with re for both active and passive net- 
works, as shown in the inset Fig. 5. At larger bending 
rigidities a c saturates to a value independent of re. Inter- 
estingly, the values of a c for active floppy networks are 
substantially lower than for the passive networks. This 
indicates that internally generated motor stress is more 
effective in network stiffening than an external stress. 

To identify the role of filament length in motor gener- 
ated stiffening we vary Q to tune (L) . Interestingly, only 
networks with relatively short filaments stiffen strongly 
(Fig. 6). Networks with longer filaments are governed 
increasingly by the stretching modes in the system. This 
is consistent with the numerical data in Fig. 5, for which 
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FIG. 6: The linear shear modulus G as a function of motor 
generated stress jm for various values of (L) at fixes bending 
rigidity k = 10 -3 . The stiffening curves for (L) < 5 show an 
approximate scaling behavior given by K ~ a, as shown by 
the dashed lines that indicate a slope of 1. The inset shows 
the critical stress for the onset of stiffening as a function of 
k for both the active (red squares) and the passive (black 
circles) systems. 



we observed that only bending dominated networks are 
capable of stiffening by motor activity. The critical stress 
for the onset of stiffening scales in the same way with (L) 
for the active networks as for the passive networks (inset 
Fig. 6), similar to what we observed for the scaling of 
cr c with k (inset Fig. 6). Taken together, these results 
provide evidence that the motor generated stiffening in 
the active networks derives from the same origin as the 
stiffening of passive networks under external shear. 

The analogy between external stress and motor gen- 
erated stress can be further explored by determining 
the effect of motor activity on the microscopic defor- 
mation field. The stiffening in passive networks has 
been attributed to a shear-induced cross-over between 
soft bending modes and stiffer stretching modes; con- 
comitant with this cross-over the deformation becomes 
increasingly affine for larger strains [29] . Our simulations 
suggest that the same basic mechanism is responsible for 
the motor generated stiffening in non-affine networks. To 
further test this picture we investigate the microscopic 
deformation field of the these networks under a small ex- 
ternal shear. We subtract the affine deformation 6r ■ of 
a cross-link i from the actual deformation ^r^ to isolate 
the non-affine contribution, 
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Consistent with prior work [24] for a passive networks 
deep in the bending dominated regime, we observe large 
non-affine deformations, as shown in Fig. 7a. In contrast, 
when motors are present the non-affine contribution to 
the deformation field is substantially reduced, as shown 
in Fig. 7b. Note, that the motors will initially generate 





FIG. 7: The non-affine component of the deformation field un- 
der an external shear for a passive (a) network (k = 10 -3 and 
Q — 1/4) and for the same network with motors (pm = 0.061, 
fo ~ 10 -2 ) deep into the stiffened regime. The greyscale of 
the arrow heads indicate the magnitude of the non-affine de- 
formation; black indicates a large magnitude (~ 0.01) and 
light grey a small magnitude (~ 0.001). The motors are shown 
as red dumbbells. 



highly non-affine deformations and large bends. These 
results show, however, that the subsequent deformation 
of this active network under a small external shear is 
considerably more affine than in the passive case. This 
provides insight into the motor induced stiffening we ob- 
served in our simulations (Figs. 5 and 6). Motor activ- 
ity pulls out the floppy bend modes, which renders the 
network deformation more affine and, thereby, induces a 
cross-over from a response governed by bending modes 
to a response governed by stretching modes. 
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III. CONCLUSION 

Here we have show that molecular motors — modeled as 
force dipoles — stiffen non-affine networks. Interestingly, 
we find that only networks that are strongly governed by 
bending modes arc capable of stiffening through motor 
activity. The internal stresses generated by the motors 
pull-out the floppy bending modes in the system, leaving 
the stiff stretching modes. In this way, motors induce a 
cross-over to a stretching dominated regime, in analogy 
to prior results on externally-stressed networks [29, 32]. 
The absence of motor-induced stiffening of our networks 
in the stretching dominated regime can be attributed to 
the purely linear force-extension behavior in our model. 
Analytical studies based on affine stretching dominated 
networks have shown that motor activity can lead to stiff- 
ening when the expected non-linear force-extension rela- 
tion is taken into account [18-20]. 

Nevertheless, within the model we consider, with 
purely linear elements, our results support the qualitative 
equivalence of external and internal stress in the nonlin- 
ear network response [4, 7, 18]. So far, this correspon- 
dence has been understood in the context of stretching- 
dominated networks, with nonlinear filaments [18-20]. 
The present work shows that this analogy is more gen- 
eral. Interestingly, however, there are some quantitative 
differences between network stiffening by external load 
vs internal motor stresses. Specifically, our results sup- 
port the idea that motor stresses can be more effective in 



generating stiffening, since they act in all directions [7]. 
By contrast, when a network is externally sheared most 
stress is focussed on a small fraction of the filaments that 
are oriented the direction of extension. Furthermore, 
there are quantitative differences in the form of the stiff- 
ening response with stress in the present model. We find 
that motor contractility leads to an increase in the shear 
modulus with motor stress ctm (Figs. 4, 5) that is approx- 
imately given by G ~ af^, where x w 1. By contrast, the 
stiffening by external shear exhibits a more complex de- 
pendence on the stress, with two distinct regimes, corre- 
sponding to x ~ 1 and x ~ 1/2. One important difference 
that sets the passive networks apart, are the geometric 
effects that arise at large external shears through the col- 
lective alignment of filament in the direction of maximum 
extension. 

The results presented here provide further insight into 
the mechanisms available for the active cellular cytoskele- 
ton to regulate the mechanical behavior of the cell. Fur- 
thermore, these principles can inspire the design of novel 
active biomemetic materials with tunable elastic proper- 
ties. 
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